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We consider A'' driven two-level atoms interacting with a structured reservoir. By dressing the 

collective operators within a semiclassical approach, we derive a master equation and a mean-field 

single particle effective Hamiltonian. This Hamiltonian describes the optical bistability phenomenon 

'_^' ' occurring in the relation between an input electromagnetic field and the effective output generated 

I I by the N atoms. The dissipative part of the master equation and the effective single-particle Hamil- 

^— V , tonian contain new terms due the reservoir structure of modes. In plotting the output field amplitude 

^vj ' and phase, for a structured reservoir, as function of the input amplitude, one verifies the bistable 

' behavior in both. We illustrate our results for two structured reservoirs: one having a Lorentzian 

f^ , shape for the distribution of modes, and the second is modeled as a photonic band-gap structure. 

< 

(N 



PACS numbers: 42.65.Pc, 42.65. Sf, 42.50.Ct 



Y^ : I. INTRODUCTION 

^ |. Optical and dynamic properties of atoms coupled to dissipative environments with a tailored density of modes have 

O ■ been a hot topic in current research |l|-[l3|. Several interesting and potentially useful effects such as: (i) suppression 

c/^ ' of spontaneous emission Ml, m, \l^ ; (ii) modifications in the resonance fluorescence and absorption spectra of strongly 

O ' driven two-level atoms [3|-[l0j , (iii) amplification without inversion 9] , and (iv) the possibility of effective control of 



c/j I atomic states [11|, have been reported by considering the interaction of atoms with structured reservoirs. In general, 
^~>. theoretical studies of dissipative dynamics in structured reservoirs has been essentially devoted to systems of one or 
i-C ' a few atoms. However, for a large number of atoms, one interesting phenomenon is the one related to the optical 
bistability (OB) |14l4l6l |. It originates from a nonlinear relation between the intensity of an input field and an effective 
output field emerging from a collection of two-level atoms. Graphically, one sees the occurrence of a S-shaped curve, 
that corresponds to the existence of two stationary stable states for the atomic system, thus allowing the use of the 
system as an optical switch [l6'|. 

Since its prediction and observation, in the 1970's, the OB |17- 19i has been the object of intense research due 



^2^ ■ its theoretical and experimental usefulness for studying nonlinear 20l-[24|. together with far from equilibrium effects 

vQ ' arising in complex systems |25l427| . The main motivation was the potential applications in the construction of optical 

• . devices ([lg| and references therein). The basic standard successful description of the OB consists in considering a 

^r ' system of homogeneously broadened two-level atoms driven by a coherent resonant field [2gl. As a matter of fact, the 

Z^ , phenomenon arises from to the interplay between (1) dissipation due to the presence of a reservoir responsible for the 

I ■ atomic decay, (2) feedback due to the mean-field many-atoms effects, and (3) pumping of the atomic sample by an 

• • I external field, occurring simultaneously. In a more specific context, in Ref. [29| the authors study the change in the 

. ^ • bistable S-shaped curve occurring in the population difference in impurity two-level atoms in a pseudophotonic band 

S^ gap background to an applied laser field. So, the structure of modes of the reservoir may influence significantly the 

^H ■ behavior of the bistable system, as we are going to further explore in this paper. 

_ Cy _■ We shall treat the problem of atomic bistability in a two-level N-atom system interacting with structured reservoirs. 

Formally we develop our calculations using the semiclassical dressed atom approach for the collective operators, 

extending the treatment considered in [3,[31 for a single atom. We derive a master equation for iiT (1 < K < N) atoms 

and for a dilute system we obtain an effective single-particle nonlinear Hamiltonian for a representative particle. The 

effective Hamiltonian and the dissipative terms in the master equation contain additional terms that are absent when 

a structureless reservoir is considered. The relation between input and output fields is obtained and, in contrast with 

the case of a structureless reservoir, we found that besides the output amplitude, also the phase presents a bistable 

behavior. So, this feature allows probing the presence of a reservoir having a non-flat structure of modes. We illustrate 

our results for: (1) a reservoir having a Lorentzian shape for the modes distribution, and (2) a reservoir displaying a 
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photonic band-gap structure. 

The article is organized as follows: In Sec. II we introduce the Hamiltonian describing the system of atoms plus 
fields. In Sec. Ill we obtain a master equation for the N atoms interacting with a structured reservoir. In Sec. IV, a 
mean-field approximation is developed and an effective single-particle Hamiltonian is obtained. In Sec. V the relation 
between input and output fields in the stationary state is obtained. In Sec VI we illustrate the bistable behavior 
considering two illustrative models for the reservoir. Finally, in Sec. VII we present a summary and conclusions. 

II. ATOMS-FIELD SYSTEM 

We consider N two-level atoms, with transition frequency ujq, interacting within the rotating wave approximation 
(RWA) with a laser, assumed as a classical electromagnetic field of frequency lul and with the electric component 
EinC^^ , having an arbitrary phase ip. Besides, the atoms interact with a reservoir at OK, made of a continuum of 
modes, which is responsible for the atomic decay. The Hamiltonian of the whole system is given by 



where 



H^Ha + Hu + Var (1) 



Ha = ^Sq + hF {e'^^^'S^ + e--''^^'S+) , (2) 

Hr = f duj {huj)b+{uj)h{uj), (3) 



Var = h duj g{uj) [b{uj)S+e'f -t- b+{uj)S-e~''^] . (4) 

The Hamiltonian ^ represents the N two-level atoms pumped by the laser field, with coupling constant F = jiEin 
(we assume all atoms having the same atomic dipole moment /i) . Admitting that the size of the atoms is much smaller 
than the laser wavelength, the two active levels of the atoms are described by the collective operators 

N N 

1=1 1=1 

and, ao{i) and cr-i-(i) are the Pauli operators for a single particle satisfying the commutation relations of the SU (2) 
algebra, [o'o{i),a±{j)] = zL2Sija±{i) and [a+{i),a-{j)] = Sijao{i). We remind that a+ = \e) {g\, cr_ = I5) (e|, 
^0 = [|e) (e| — \g) {g\] /2, where |e) and \g) refer to the higher and lower energy levels, respectively, while hajo is energy 
difference between the levels. The Hamiltonian ^ represents the reservoir modes, where the operator b{Lu) (&+(cli)) 
annihilates (creates) a quantum of frequency w, and both satisfy the bosonic commutation relations, [6(0;), 6+ (cj')] = 
(5 (w — w'). Finally, Hamiltonian ^ corresponds to the coupling between the reservoir modes and atoms, and g{oj) is 
the coupling parameter, assumed to be frequency dependent, that characterizes a structured reservoir. Furthermore, 
we assume that the atomic system is quite diluted such that we disregard the direct interaction between the atoms, 
so they will correlate and feel each other indirectly, as an effect of their coupling with the reservoir modes. We also 
consider an undepleted laser field so its dynamics is not taken into account. 

In a referential frame rotating at frequency ujl and in the interaction picture, with respect to the reservoir modes, 
the Hamiltonian ([T]) becomes 

H = Hos + VAR{t), (6) 

where 

Hos^SSo + F{S-+S+), (7) 

and S — (wq — (^l)/2 is the detuning frequency between atoms and laser field; and 

VAR{t) = j diu g{ijj) \b{ijj)S+e'^"^-'^^'+"^ + h.c] . (8) 



is the interaction between atoms and the reservoir. 

In order to derive a dynamical equation for the atomic system, usually a master equation, we follow an approach 
considered in references [7|, |8| . Instead of considering the atomic decay process independently of the driving field, one 



assumes a coupled atom-driving field decay, i. e., a generic atom is dressed by the driving field and coupled to the 
reservoir modes. Within this approach we define the semiclassical dressed collective operators 



1 



5o = X (-5^0 + F{S- + S+)) 



S+^^ {{S - A)S^ - FSo + {5 + A)S+)) ., 



(9a) 
(9b) 
(9c) 



that satisfy the very same commutation relation of the SU (2) algebra [S'o,S'±] = ±2S+, [S+,S-] ~ Sq. For later 
convenience we write Eqs. (|9a| -(|9c |) in short as 



Sj == ^ CjiSi, j = -1,0,1 



(10) 



i=-i 



where Cji are the entries of the matrix 



•iF,S,A) 



2A 



S + A -F S-A 

2F 25 2F 
6- A -F 5 + A 



(11) 



and A = \/5'^ + F"^. Definition (|10p is invertible, so the operators S± are related to the dressed ones through the 
relation 






(12) 



As C^'^ (F) = C (— -f), the entries of the matrix C {F, S, A) are Ckj (F) = Ckj {—F), and the dependence on the other 
parameters remains the same. 

Defining the single index coefhcients Cj — cij — C-i-j (where Cj=±i — {5 + jA) / (2A) ,Co — F/2A), the 
interaction Hamiltonian ([8]) can be written as 



VABit) = J2^n ^^ 5(^) [Bj {t,u) Sj + b] {t,io) 5_ 



(13) 



where Bj (t) = b{Lu)e 



-i{uj—ujL—2A.jt)t-\-iip 



III. iV-ATOM MASTER EQUATION 



In order to describe the dynamical evolution of an A'^-atom system state we begin with the general non-Markovian 
master equation in the interaction picture [SOi] 



^^^f^ ^-Jd^ Tr7^ [VARit), [F^fl(i - r), p,,^(i - r)p^]] 



(14) 



for the evolution of the density operator /O/, at (t) of N atoms coupled to an unperturbed reservoir in thermal equilibrium. 
As usual, Eq. ([T4| is obtained by tracing over the reservoir degrees of freedom, and pR is the reservoir density operator. 
The index / stands for interaction picture. Besides, Eq. (J14p was obtained in the weak coupling approximation, by 
assuming that the atomic system plus reservoir density operator factorizes as ps+R{t) = ps{t)pR, at any time. 
By inserting the interaction term given by Eq. (J13p in the master equation (|14l) we obtain 



dpi,N{t) 
dt 



1 1 /'OO pC 

Y, Y, C,Cy \ du^gicj) / 



duj' g{Lo') I dr 



xTr 



n 



[s, (i, w) Sj + b] {t, w) S^j , \Bj, (t - r, u) S,, + b], {t - t, w) S.j, , p/,jv(i - r)p^ 



(15) 



We note that the density operator of the iV-atom system at time t depends on the density operator at the previous 
time t — T. So, at this point, we invoke the Markov assumption by replacing pi^N{t — t) in Eq. (|15p by pi^N{t) and 
extend the upper hmit in the integral to infinity. As the reservoir is assumed in a vacuum state (OK)^ within that 
approximation the trace operations result in 



TrTz{b{^)b\u>')pR} = 5{ll>-uj), 
Tnz{b\Lo)b{uj')pii} - 0, 



(16) 
(17) 



and Eq. (ITSI) becomes 



dt 

where 



E ^^^T 



JiA{j+ji)t 



Cj, (c^L + 2jA)([^,,5_,vp/,jv(t) 



^2iA{j-]i)t 



S-j,Pl,N{t)Sj 



'2iA{]-]i)t 



/"OO 

^j {ujL + 2jA) = / duj (giiojf 
Jo 



TrS{Lu ~UJL- 2jA) ± iV 



1 



cj — ujl — 2j A 



(18) 



(19) 



and V stands for principal part. The first term at right hand side of Eq. (|T9|) is responsible by the atomic decay rate 
while the second term gives a shift in the atomic frequencies. Usually, these shifts are incorporated in the atomic 
frequencies, but since our aim is to study the decay rates, in this work it will be neglected and we write 



^,^ng\u;L + 2jA) = ^ 



(20) 



(j = —1,0, +1), where we introduced Tj as the decay rates that depend on the driving field and V is the volume of 
the atomic cell that contains the N atoms. 

Back to the Schrodinger representation, the master equation for an A'^-atom system becomes. 



dt 
fo 



= -i[HQs,PN{t)] - 7rr7{[S+,S-PN] +H.C.} 



{[S+,SqPn] +h.c.} 



2V 



{[S+,S+pn] +H.C.}, 



where H.c. means Hermitian conjugate, and the effective decay rates are given by 



To = 



1 



4A2 



-F2 2^2 _p2 

-{A + 5)F 2SF {A-S)F 
(A + 5)2 2F^ (A -(5)2 




(21) 



(22) 



So, in a structured reservoir the A^-atom system contains decay rate parameters Tj that have a dependence on the 
frequency and the intensity of the laser field. In contrast, in a non structured reservoir, g{io) = go, the decay rates 
are either constant or zero, r+ = 27rg2 = F, Fq = r_ = 0, and the master equation (|2T|) reduces to the well know 
form [30i] 



dpNJt) 
dt 



'i[H(iS,PN{t)] 



2V 



{[S+,S^Pn]+H.c.}. 



(23) 



The extra terms, third and fourth, appearing in Eq. (j2ip look like the dissipative terms that appear due to a squezeed 
vacuum reservoir [32| and where obtained in the context of a single atom in [8] . 



IV. 



MEAN-FIELD APPROXIMATION AND EFFECTIVE HAMILTONIAN 



Now we will treat the A^-atom system as a quantum BBGKY hierarchy of equations similar to that of classical 
kinetic theory [3a, [Sg. We consider a subsystem constituted of K atoms, with K < N. The equation of motion for 



the K atoms density operator pK is obtained by calculating the trace over the remaining K + l.K + 2, ...,N atoms 
degrees of freedom in Eq. (PT|) , so getting 



dPK „ { dpN 

— J rK+i N 



dt 



dt 



-iTrK+i,...,N {5[So, Pn] + F[S+,pn] + F[S-,pn]) 



2V 



{TiK+i,...,N {[S+, S-Pn]) + TrK+i,...,N {[pnS+,S-])} 



- -Z7T {TrK+i,...,N i[S+, SoPn]) + TrK+i,...M {[pnSi^, S-])} 

- -^ {TrK+i,...,Af {[S+, S+pn]) + TrK+i,...,N {[pnS-,SJ\)} . 
Using the sums ([Sj in terms of the microscopic operators Eq. ([24|) becomes 

K 



dp 



K 



dt 



= -*^ {[5ao(i) + Fe-'''Pa+{i)+Fe''Pa-(i),pK]] 



i=l 



K 



[N ^).^ J Lt[^^(,)^TrK+i(a_(if + IW+1)] 



V 



(24) 



+ ^e-na+(z),TrK+i (fTo(if + iW+i)] 

+ £_e-2^^[a+(i), Tr^+i {^+{K + 1)pk+i)] + H. 



K 



-E 



2,2' — ! 



^[a+(z),a_(z')pA-] + ^e-^[a+(^),ao(^')pK] 



+ — e-^'^^[a+{i),a+{i')pK]\+H.c. 



(25) 



For a single representative atom of the system {K = 1), the equation of motion for pi will depend on the two atoms 
density operator p2, and so on for the whole hierarchy, meaning that the equation of motion for pK will depend on 
the state pk+i- However, for a dilute system the higher-order atomic correlations may be disregarded, and this is 
achieved when one factorizes p2 as pi® pi, so a generic atom is assumed to move in a mean- field produced by all the 
other atoms, which is a kind of Hartree approximation. Implementing this approximation and dropping the subscript 
in pi, the equation (|25p reduces to 

— = - i[Hef [p],p\ ~ ^ ([cr+,CT_p] + H.c.) 
- ^{e-'^[a+,a,p]+H.c) 



— (e-^^^[c7+,a+p]+i/., 



where the single particle effective Hamiltonian is 



(26) 



Hef [P] = 5(Tq 



pEi, 



AN-i) 

V 



{a+)e 



{<^o) 



H.c. 



(27) 



which contains nonlinear terms corresponding to a mean-field that is due to the remaining iV — 1 atoms. 

Prom the second term in the brackets in the Hamiltonian (l?71) we see that effectively a single generic atom is excited 
by the input field amplitude £"„ plus an extra polarization field density e* /(i), where 



epoiit) = »^^^^e'^ (f_(a_)e'^ + f o(ao) + f +(a+)e-^^) 



(28) 



originated from the other {N — 1) atoms that produce a mean-field effect, and is proportional to the uniform atomic 
density in the cell N/V, for iV ^ 1. When the function g^(w) is assumed being frequency independent (white noise 
limit) the polarization field reduces to the simple expression 

epoi{t) = t^-^^^-^T{<j+}. (29) 

The equations of motion for the atomic operators mean values, derived from the master equation (j26p . are 

-^(ao) - 2i^i {eout{t){a^) - eo„t*(i)(a-)*) + fo ((a+)e-"^ + (a_)e'^) - f+(l + (,so)), (30) 



dt 



-(a_) = -2z5(a_) +*Meo«t*(t)(cTo) + ^ (f_(a+)e-*'^ -f+(a_)e'^) + -^e'^^ (31) 

and (cr+) = {(T-Y ■ The total effective output field transmitted from the sample is defined as 
f-out{t) = E^ne'-'^ + epoi{t) 

= il;„e*'^+^^^^^f+(a+)+^i^^-^e''^(f_(a_)e^^ + f^ (32) 

When g{uj) — go, in Eq. ([32]), the only remaining term is r+ — > F (proportional to cr+), a constant, while the other 
terms vanish. The mean- field extra terms, proportional to (cr_) and (uo), are due to the structured reservoir. The 
Eq. ([5^ . for the output field tout{i), is more inclusive than the others deduced in Refs. [2^ [2^ [SJ], because the 
mean-field approximation contributes with additional terms that are sensible to the mode-structured reservoir. 

V. STATIONARY SOLUTION AND THE INPUT-OUTPUT FIELDS RELATION 

When the solutions of Eqs. ([50]) and pil) are inserted in Eq. ([5^ we get the output field (total field) as a function 
of the input one. Here we are interested to study the influence of the structured reservoir on the bistable steady state 
output field amplitude as a function of input field £'i„. The stationary solutions of Eqs. (I30p and pip are obtained 
by setting d{ao)/dt — d{a^)/dt = 0, resulting in 

, , _ f+(16(52+f2.-f^)-2f2(f++f_)-2fo[i(-4i(5 + f++f_)/xe,,e-*'^ + c.c.] 

(Cq/ss — ' — ~ ~ ~ ~ ~ ~ ~ ~ o 

r+(16(52 +Tl- Tl) + 2nTQ[i{AiS + T++ r_)es.e-'¥' + c.c] - Vr_ (e2^e-2iv -|- c.c.) + 8^i^T+ |e,,r 

(33) 
and 



-'^(^+)s. 



{-f of+ (aiS + f+ + f_) + 2i[tl - T+iAiS + ^+)]^lesse~ 



■l^ 



r+(1652 + r2_ - Tl) + 2fiTo[i{4iS + r+ + r_)esse-'-fi + c.c] - Vr_ (e2,e-2''^ + c.c.) + 8iJ.^r+ \ess\ 

(34) 
where e^s is the output field at the stationary state. 

In particular, for a non-structured (ns) reservoir, g{uj) — go, Eqs. ((33)l and (p4)) simplify to the well known results 

M 

iao)^^ = (^^^^ + ^^) „ e-^(a,)(r>=^^-^^ + -^)^----". (35) 

16(52 + r2+8Ai2|e,,|' "' 16(52 + r2 + 8^2 |g^^|2 

for a single atom pumped by an external field Ess = EinC^'^ + i {{N — l)/2fiV) T{(t+)ss, which is due to a collective 
effect produced by the mean-field. For a structured reservoir the stationary solution changes in an essential way: 
besides of the terms proportional to \ess\ , others terms, proportional to ess and £^„, appear additionally. These terms 



show some similarity to those produced by the decay in a squeezed vacuum [3l|-|3J|, without mean-field effects, as 
reported in |8| in the case of a single atom. 



Writing ess = EssC^^ (with the exphcit introduction of a phase 9) in the stationary state of Eq. (|32|) . we obtain the 
following relation between the input and output fields 



Et„ = \EssCos<i> 



JN-l) 
2Vfi 



f+-f_)lm(e-'^(a+)ss)) 



i^..sin$-(^(r, + r_ 



Re e-'^(CT4 



ro(cro) 



(36) 



where {a-\-)ss and (cro)ss, from Eqs. (p3| and (p4l) . depend on the phase difference ^ = 6 — ip. We note that the 
nonlinear dependence of 9 on Ein is a manifestation of the intrinsic frequency distribution of the reservoir. Thus, 
in order to determine graphically the relations between the real and imaginary parts of ess versus Ei„, or Ess and 
$ versus E^, for each input value Ein, the two output dependent variables must satisfy Eq. ([55]) for the unique 
independent variable Ein- In the white nois e ap proximation we have r_ = Fq = 0, F^ = F, so Eq. p6p simplifies to 
the already known phase independent form [l6[, 



Ei, 



Es 



(N-l) 



1/(1 + 16 (<52/F2) + 8fi^El/T^) 



4:{N -i)s/r 



y (1 + 16 {syr^) + Sfi^Ei/r^) 



1/2 



(37) 



which displays the bistability phenomenon when output versus input fields are plotted. The first term in the braces of 
Eq. p7p corresponds to the absorptive regime, when the atoms are driven near resonance, while the second term stands 
for the dispersive regime, when the atoms are driven far from resonance (S/T ^ 1) and nonlinear refractive effects 
dominate [16|. Very characteristically, in this approximation (structureless reservoir), there is no phase dependence 
induced by the atoms. Therefore the presence of a phase shift indicates the existence of a structured reservoir (colored 
noise). 

VI. STRUCTURED RESERVOIRS: OUTPUT FIELD AMPLITUDE AND PHASE BISTABILITY 

In order to get a direct insight into the problem of bistability with a structured reservoir we consider the case when 
frequencies of the laser field and atomic transition are resonant (S = 0). Furthermore, calling Re(eout) = Ex and 
Im(eoMt) = iy, Eq. (I55| becomes 



(^o; 



(f+ + f_) [f+ (f+ - f_) - 2f 2 + 4^Itoe, 



F. F, -r_ 



VFoEy + 8Ai2e2 



8^2 r+-r_ e; 



and the real and imaginary components of ((T+)^^ are, 

(fo - 2^iey^ f+ (f+ + f_) + Sfi^tot 



Re (aA 



T. T 



ifiToey + Sfi^el 



8^,^(i 



r\e 



(38) 



(39) 



Im(o'+).s = 



F . F 



2F2 + i^iVoey 



M^x 



f+ + f_) [f+ (f + - f_) - 4^foe, + Sfi^l 



8/i2 F+ - F_ e 



(40) 



Inserting Eqs. pSMOp into Eq. ([5^ . and equating the real and imaginary parts, we get the following set of 
equations. 



£;. 



tn — ^x 



2Vn 

JN-l) 
2Vfi 



F+-F_ Im((a+),,) 



F++F_ 



Rc{{a+)ss)+fQ{(To), 



(41) 
(42) 



for the components ex and e^, and whose solutions determine Ess — ^/e^ + ^y and $ — arctan(ej,/ea;). 



By assuming, in particular, that the function g'^{uj) is symmetric around the atomic transition frequency - an even 



function with respect to ojq, g {ojq 



g {ojo + uj) -, the decay rates r_|_ and T_ become equal, so the effective 



decay rates in Eq. ([22]) reduce to f ± = ^ (Fq ± r+) and f o = 0. In this case the output field does not depend on 
(ctq), from Eq. P^ we have ey = 0, and the unique physical solution occurs for $ = as in the case of a structureless 
reservoir. Nevertheless, the decay rates depend on the input field, so the relation between input and output field 
amplitudes simplifies to 



Ej, 



N-i r+ (^„o 



V 



To 



ro+r+(B,„) 



(43) 



where we wrote explicitly the dependence on £*,-„,. In the case of a structureless reservoir (white noise), Fq = r+, Eq. 
(|^5)l reduces to the well known equation [1^, ^^ 



Ej, 



E. 



(iV-1) 



V {l + ^^i^ElJT^) 



(44) 



whose quite simple input-output amplitude nonlinearity is due to the term proportional to iV— 1. In order to compare 
with previous works [22.], using the authors' notation, we write Eq. (|43p as 



y = x- 



'iC{y)x 
D{y)+x'^ 



(45) 



where we have defined the input and output fields as y — ySfiEin/To and x = y8fj,ET/^o, respectively, C{y) — 
((7V-l)/y)(Fo + F+)F+/4F2 andi:>(y) = (Fo + F+)/2Fo. For a structureless vacuum we have C(2/) = C= {N-l)/V, 
D{y) — 1, so Eq. (|i5|) reduces to 



y = x + 



2Cx 

1 + a;2 



(46) 



that is identical to the well known result obtained in 28]. Equation (|46)) displays a bistable behavior for C > 4, and 
for C ^ 1, the range for input field, allowing three solutions, lies in the interval \/8C < y < C. Now, regarding Eq. 
()45[) . the dependence of x on y, is obtained by inverting the expression and solving the cubic equation, which, for 
a structured environment, has its coefficients depending nonlinearly on the input field amplitude. So, the range of 
values for a bistable solution is determined by the density of particles and by the proper input field. 

For a more general physical situation, when g{uj) is not symmetric around the atomic frequency loq, all terms in 
Eq. (p2|) contribute to the output field. In this case Tq ^ and an imaginary component ey arises in the output field. 
Thus, in contrast with the case of a symmetric g{uj) (or just being fiat), the asymmetric structure of modes induces 
a relative phase in the output field, and both, amplitude and phase, display a bistable behavior. We illustrate below 
the input-output field dependence for two reservoirs models. 



A. Example 1: weighted Lorentzian shape 



We assume a structured reservoir having a weighted Lorentzian shape for the frequency distribution [y, [l/ 



2/ N ^0 

Ztt 



/3+(l-/3) 



7 



(a; - ujr)^ + 7^ 



(47) 



with ojr is a characteristic frequency of the reservoir and the parameter < /3 < 1. In Eq. (H71) . the first term in the 
brackets represents a background vacuum (white noise) whereas the second term represents the structured vacuum 
(colored noise), assumed to have a Lorentzian shape of width 7; the parameter (3 interpolates between the two limiting 
cases respectively. The effective decay rates are 
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where rj = tUr — ujq is the detuning between atomic frequency and the reservoir characteristic frequency. Worth to 
note that while T_ and r_|_ are even functions regarding the change 77 — > —77, Fq is an odd function. The existence 
of an input-output phase difference depends on a nonzero Fq, and when one changes the sign of rj, also that phase 
changes sign. 

In Figs. [T]we plotted the output versus input amplitudes for three values of /?, a detuning j] — 0, and N = 50. In 
this case ej, = and e^ = Ess, so the phase difference is $ = 0. In both cases, /? = 0.5 and /3 — 0.0 we observe a 
variation in the distances between the switching points, P and Q, that indicate the location of the lower and upper 
branches of the S-shaped curve. Comparing to the structureless reservoir, /? = 1.0 (solid line), there is a reduction in 
the range of values of £"„ where the bistability occurs. We also observe that the deviations from the white noise curve 
(solid line) is more pronounced at the switching points (P) from lower to upper branches. As the effective decay rate 
diminishes with the increase of the input field, less energy is transfered from the atomic sample to the reservoir, thus 
the energy goes through the sample carried by the output field. We now consider the detuning 77 = 5.0 and draw in 
Fig. [5^ and [^b the output phase and amplitude, respectively, showing bistability in both. In Fig [5^ we observe that 
the phase difference is not null because of the emergence of the component e^ . The variation of $ is more pronounced 
the more the mode distribution (/3 = 0.5 and /3 = 0.0) deviates from the white noise {(3 = 1.0). The phase $ goes to 
zero for large values of the input amplitudes because the effective decay rate Fq goes to zero. Regarding the /3 < 1 
cases, we note that amplitude and phase have the same switching points, and the P - Q distance is reduced. By 
admitting a negative detuning 77 = —5.0 the input-output amplitude relation does not change, however, as can be 
seen in FigOthe bistable behavior of the phase changes by a sign inversion. 

B. Example 2: photonic band gap 

Here we adopt a simple reservoir structure used to analyze the resonance fluorescence phenomenon in a photonic 
band gap [10] , where it is assumed that there is a discontinuity at specific frequencies of the photonic density of modes 
17^(0;), although it is constant over spectral regions in the dressed atomic frequencies; so 
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In this case, at each dressed frequency, i.e. Ein, ^ 0, the decay rates are F_(ajo — 2^Ein/fi) = Fi and Fq (ojq) = 
F+ {luq + 2^iEin/h) = F2, and the effective decay rates become fo = f _ = (F2 - Fi) /4 and f+ = (Fi + 3F2) /4. 

In Figs. |4^ and Hb we plotted the output field phase and amplitude, respectively, as a function of the input 
amplitude, for some ratios of F2/F1. We set N = 50, and all parameters are dimensionless. In Fig 0^, as in the 
previous example, we observe a change in the amplitude input range for bistability. Also, the switching point from 
lower to upper branch is more sensitive to the ratio F2/F1. Regarding the phase, from Fig. |41d we observe that, 
depending on the ratio F2/F1, at the region of bistability the phase difference may change sign and also can display a 
more complex behavior (a loop) than an S-shaped form. For large values of the input amplitude, from Eqs. p8ll42p . 
we get ex = E^, ey attains, asymptotically, the constant value {N — 1)F2(F2 — Fi)/4F^(Fi + F2), and the phase 
$ goes to zero with a sign that depends on the difference F2 — Fi. For larger values of Ein the phase difference $ 
becomes proportional to the jump (discontinuity) in the structure of modes. 

VII. SUMMARY AND CONCLUSIONS 

We presented a study of optical bistability using the mean-field approximation by considering a system of N two- 
level atoms interacting with structured reservoirs. Our approach consisted in dressing the atoms collective operators 
with the classical input field and coupling them to the structured reservoir. The dynamical system is described by 
a master equation which contains extra terms (compared to that obtained from a structureless reservoir) resembling 
those present in the master equation, derived under the influence of a squeezed reservoir. The master equation contains 
effective decay rates, associated to each dressed frequency, that depend on the frequency and intensity of the input 
field. Adopting the mean-field approximation, we deduced a single particle effective Hamiltonian containing extra 
nonlinear terms which are absent in the case of a structureless reservoir. In the stationary regime of the atomic and 
laser field system, we analyzed the relation between the input and output fields, observing that they are related in a 
nonlinear form and present typical bistable behavior. For a structured reservoir the S-shaped curve is present in the 
amplitude, but not necessarily in the phase. This theoretical result compares with a similar treatment given in [8| 
where the authors analyzed the resonance fluorescence and absorption spectra of a single atom and in which a similar 



10 

phase dependence in the system dynamics occurs. However, both the resonance fluorescence and the absorption 
spectra have no phase dependence. Our results indicate that the presence of the induced phase shift in the output 
field of an N-atom system could be an interesting probe about the nature of the reservoir. We have considered in 
details the case of resonance between atoms and the input field frequencies and verified that the output phase shift 
appears for reservoirs having an asymmetric structure of modes. We presented two illustrative examples of reservoirs: 
(1) a mixing of white noise and Lorcntzian shaped frequency distribution and (2) a simplified photonic band-gap 
distribution of the reservoir modes. We noted that the output field bistable behavior is as sensible in the phase 
difference, or acquired phase, as for the amplitude. That characteristic can be explored by using a Mach-Zehnder 
interferometer, with a phase shifter in one of the arms, where the output field phase can be determined by measuring 
the difference of the pulses intensities at the interferometer output ports. By slightly changing the input field Ei^ 
at the transition values in the instability region, it should be noted, as a response, the occurrence of sudden and 
discontinuous changes in the difference between the photocurrents at the two exit ports of the interferometer. On the 
other hand, by engineering reservoirs and atomic samples one could use the output amplitude and phase to design 
optical control devices. Besides, we will address, in a future work, the influence of structured reservoirs on bistability, 
in the dispersive regime (S ^ 0). 



Acknow^ledgments 

GAP, ACY, and SSM, acknowledge financial support from FAPESP (Sao Paulo, SP, Brazfl), through the INCT-IQ, 
and CNPq (DF, Brazfl). 



[1] 


M 


[2] 


T. 


[3] 


J. 


[4] 


J. 


[5] 


J. 


[6] 


G. 


[7] 


C. 


[8] 


A. 


[9] 


M 


[10] 


M 


[11] 


G. 


[12] 


P. 


[13] 


D. 


[14] 


H. 


[15] 


L. 


[16] 


A. 


[17] 


A. 


[18] 


S. 


[19] 


H. 


[20] 


H. 


[21] 


K. 


[22] 


H. 


[23] 


M 


[24] 


H. 


[25] 


C. 


[26] 


F. 


[27] 


F. 


[28] 


R. 


[29] 


S. 


[30] 


D. 


[31] 


S. 


[32] 


c. 


[33] 


J. 


[34] 


L. 



Lewenstein, T.W. Mossberg, and R.J. Glauber, Phys. Rev. Lett. 59, 775 (1987). 

W. Mossberg and M. Lewenstein, in Advances in Atomic, Molecular and Optical Physics, edited by P.R. Berman, 
Academic, New York,1994, Suppl. 2, p. 171. 

R. Brinati, S. S. Mizrahi, and G. A. Prataviera, Phys. Rev. A 50, 3304 (1994). 
R. Brinati, S. S. Mizrahi, and G. A. Prataviera, Phys. Rev. A 52, 2804 (1995). 
R. Brinati, S. S. Mizrahi, and G. A. Prataviera, Phys. Rev. A 56, 322 (1997). 

A. Prataviera, S. S. Mizrahi, V. V. Dodonov, and J.R. Brinati, Phys. rev. A 60, 4045 (1999). 

H. Keitel, P. L. Knight, L. M. Narducci, and M. O. scuUy, Opt. Commun. 118,143 (1995). 

Kowalewska-Kudlaszyk and R. Tanas, J. of Mod. Opt. 48, 347 (2001). 
, Erhard and C. H. Keitel, Opt. Commun. 179, 517 (2000). 
, Florescu and S. John, Phys. Rev. A 69, 053810 (2004). 

Harel, A. G. Kofman, A. Kozhekin, and G. Kurizki, Opt. Express 2, 355 (1998). 

Lambropoulos, G. M Nikolopoulos, T. R. Nielsen, and S. Bay, Rep. Prog. Phys. 63, 455 (2000). 

G. Angelakis, P. L. Knight, and E. Paspalakis, Contemporary Physics 45, 303 (2007). 

M. Gibbs, Optical BistabtUty: Controling light with light, (academic press, New York, 1985). 

A. Lugiato, progress in Optics, edited by E. Wolf (North Holland Amsterdam, 1984), Vol. XXI, p. 69. 

Joshi and Min Xiao, J. Mod. Opt. 57, 1196 (2010). 

Szoke, V. Daneu, J. Goldhar, and N. A. Kurnit, Appl. Phys. Lett. 15, 376 (1969). 
L. McCall, Phys. Rev. A 9, 1515 (1974); 

M. Gibbs, S. L. McCall, and T. N. C. Vekatesan, Phys. Rev. Lett. 36, 1135 (1976). 

M. Gibbs, F. A. Hopf, D. L. Kaplan, and R. L. Shoemaker, Phys. Rev. Lett. 46, 474 (1981). 

Ikeda, H. Daido, and O. Akimoto, Phys. Rev. Lett. 45, 709 (1980). 

Nakatsuka, S. Asaka, H. Itoh, K. Ikeda, and M. Matsuoka, Phys. Rev. Lett. 50, 109 (1983). 
, Taki, Phys. Rev. E 56, 6033 (1997). 

A. Babu and H. Wanare Phys. Rev. A 83, 033819 (2011). 

M. Bowden and C. C. Sung, Phys. Rev. A 19, 2392 (1979). 

On and Z. Qin, Opt. Commun. 65, 455 (1988). 

On, Phys. Rev. A 41, 3021 (1990). 

Bonifacio and L. A. Lugiato, Opt. Commun. 19, 172 (1976). 
John and T. Quang, Phys. Rev. A 54, 4479 (1996). 

F. Walls and G. J. Milburn, Quantum Optics, (Springer, Berlin, Heidelberg, 1994). 
F. Haas and M. Sargent III, Opt. Commun. 5, 366 (1990). 

W. Gardiner, Phys. Rev. Lett. 56, 1917 (1986). 
Bergou and D. Zhao, Phys. Rev. A 52, 1550 (1995). 

P. Maia, G. A. Prataviera, and S. S. Mizrahi, Phys. Rev. A 69, 053802 (2004). 



11 



[35] R. L. LibofT, Kinetic Theory, (Prentice Hall, New Jersey, 1990). 
[36] D. A. McQuarrie, Statistical Mechanics, (Harper, New York, 1976). 



12 



10 



LU 






8 - 



6 - 



4 - 



2 - 







P = 0.0 

(3 = 0.5 

(3 = 1.0 

r| = 




Hi ::^ 







8 



10 



in 



FIG. 1: Output versus input field amplitudes for rj = 0, N = 50, 7/ro = 20 and different values of /?. The solid line corresponds 
to an structureless reservoir, 13 = 1. The parameters are dimensionless. 
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FIG. 2: Output (a) phase difTerence $, in radians, and (b) amplitude Ess, versus input amplitude for rj — 5, N — 50, 7/ro = 20 
and different values of /3. The solid line corresponds to a structureless reservoir, P — 1. The parameters are dimensionless. 
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FIG. 3: The same as in Fig. [2] but for 77 = -5. 
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FIG. 4: Output (a) phase difference $, in radians, and (b) amplitude Ess, versus input field amplitudes for A*' = 50, and 
different values of r2/ri. The solid line corresponds to a structureless reservoir, r2/ri = 1. The parameters are dimensionless. 



